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Lattice Landau gauge and other related lattice gauge fixing schemes are known 

to violate spectral positivity. The most direct sign of the violation is the rise of the 

<^ ' effective mass as a function of distance. The origin of this phenomenon lies in the 

o ; 

Cn ■ quenched character of the auxiliary field g used to implement lattice gauge fixing, 

^ . and is similar to quenched QCD in this respect. This is best studied using the PJLZ 

^-5 ■ 
Q>^ \ formalism, leading to a class of covariant gauges similar to the one-parameter class 

of covariant gauges commonly used in continuum gauge theories. Soluble models 

>\ 
'nI" . are used to illustrate the origin of the violation of spectral positivity. The phase 

O \ diagram of the lattice theory, as a function of the gauge coupling (3 and the gauge- 



o 
o 



fixing parameter a, is similar to that of the unquenched theory, a Higgs model of a 
type first studied by Fradkin and Shenker. The gluon propagator is interpreted as 
yielding bound states in the confined phase, and a mixture of fundamental particles 
Q \ in the Higgs phase, but lattice simulation shows the two phases are connected. Gauge 



^ 



field propagators from the simulation of an SU{2) lattice gauge theory on a 20^ lattice 



^ . are well described by a quenched mass-mixing model. The mass of the lightest state, 

H ; 

5t , which we interpret as the gluon mass, appears to be independent of a for sufficiently 

large a. 



I. INTRODUCTION 

Although many interesting properties of a lattice gauge theory can be determined without 
gauge fixing, there are several reasons why it is needed. At a fundamental level, gauge fixing 
is necessary to make the connection between continuum and lattice gauge fields. Gauge 
fixing has also been a key technique in lattice studies of confinement as well yj, and may 
help differentiate between different models of confinement. For example, continuum theories 
of the origin of confinement often make predictions about the gauge field propagator; see 



Ref. [2] for a review. Finally, gauge fixing may be needed to determine important properties 

the q 
the finite-temperature gluon propagator 



of the quark-gluon plasma phase of QCD, such as screening masses, which are contained in 
Techniques for lattice gauge fixing have been known for some time p|- It has been clear 



from the beginning that non-Abelian lattice gauge field propagators show a violation of spec- 
tral positivity. In a theory respecting spectral positivity, the Kallen-Lehmann representation 
for a typical wall-to-wall two-point function has the form 

G{t) = dn? p (m^) e~"^l*l (1) 

Jo 2t7i 

where p is everywhere non-negative. Wall sources set all momenta transverse to the direction 
of propagation to zero. The effective mass defined by 

^e//(t) = -|lnG(t) (2) 

can be thought of as an average mass (m(t)), where the average is taken over 
p (jn?) exp(— m \t\)/2m. It follows that 

|me//(t) = (m(t))^-(m^(t)><0. (3) 

The positivity of p implies that the effective mass is a non- increasing function of t which 
goes to the lightest mass propagated in the limit t — * oo. 

Gauge invariant operators which couple to glueball states show this behavior, but the 
gluon propagator does not. Covariant gauge gluon propagators have an effective mass in- 
creasing with distance. This is not too surprising: we know from perturbation theory that 
covariant gauges contain states of negative norm. However, that knowledge has neither 
explained the form of the lattice gluon propagator nor aided in the interpretation of the 
mass parameters measured from it. In fact, the particular form of spectral positivity viola- 
tion observed for non-Abelian models is not observed in the U{1) case 8|, which also has 
negative-norm states in covariant gauges. 

In lattice simulations, gauge fixing has typically involved choosing a particular config- 
uration on each gauge orbit. A brief review of this approach is given in Ref. |2]. In the 
continuum, on the other hand, gauge fixing most often includes a parameter that causes the 
functional integral to peak around a particular configuration on the gauge orbit. The PJLZ 
formalism, a comparable formalism for lattice gauge fixing, first appeared in Refs. |9|, |lO| . 



The strong-coupling expansion was developed in Ref. [ll| and the gluon propagator was 
studied in Ref. [l^]- This formalism was used in a discussion of Abelian projection in lattice 
theories jl^. There it was shown that Abelian projection without gauge fixing leads to an 
equality of the asymptotic string tensions in the underlying non- Abelian theory and the pro- 
jected Abelian theory. Furthermore, the string tensions were proven to be equal with gauge 
fixing, provided the gauge fixing procedure respects spectral positivity. The failure of spec- 
tral positivity demonstrated in Ref. [1^ explains the string tension discrepancies noted in 
lattice studies of Abelian projection. A lattice study of the phase structure of the projected 
theory as a function of the gauge-fixing parameter was carried out in Ref. |l5l |. 

Our aim in this work is to demonstrate that spectral positivity is violated in lattice gauge 
fixing by mechanisms related to those encountered in quenched QCD. Specifically, we argue 
that the gluon propagator has a double pole structure similar to that of the r/'. We show in 
Sec. |n]that the PJLZ formalism makes clear that lattice gauge fixing is a form of quenching, 
with the gauge transformations acting as quenched fields. In Sec. IllH we develop a general 
formalism for quenched fields, and apply that formalism to gauge fixing. We show that in 
the case of a U{1) gauge theory, spectral positivity is maintained. Sec. IIVI examines two 
useful simple models for quenched fields based on mass mixing between a quenched and an 
unquenched field. The first model involves scalars; the second involves vector fields, and will 
be used to fit the results of our lattice simulations. In Sec. El we discuss the phase diagram 
of the gauge-fixed model, using the interpretation of gauge fixing as a quenched Higgs model. 
As first discussed by Fradkin and Shenker, the nominal Higgs and confining phases of the 
unquenched model are connected, and this carries over to the quenched version. We propose 
interpretations for the lattice gauge field propagator in both phases. Section |VT1 discusses 
the simulation results for SU{2) lattice propagators. A comparison of our results with some 
of the other proposed forms for the lattice gauge field propagator is performed in Sec. IVIH 
and Sec. IVIIll gives our conclusions. 

II. LATTICE GAUGE FIXING AS QUENCHING 



The standard approach to lattice gauge fixing is a two-step process |2|. The gauge fields 
U^ (x) are associated with links of the lattice, and take on values in a compact Lie group G. 
An ensemble of lattice gauge field configurations is generated using standard Monte Carlo 



methods. This ensemble of G-field configurations is generated by a functional integral 



w 



[dU] e^^t^l . (4) 



Sjj is a gauge-invariant action for the gauge fields, e.g., the Wilson action for SU{N) gauge 
fields: 

plaq 

where Upiaq is a plaquette variable composed from link variables, and the sum is over all 
plaquettes of the lattice. The gauge action Su is invariant under gauge transformations of 
the form 

Uf, (x) ^ g^ {x + /i) Uf, (x) g (x) . (6) 

The expectation value of any observable O is given formally by 

(O) = ^ I [dU] e'-^''^ O , (7) 

but in simulations is evaluated by an average over an ensemble of field configurations: 

{0) = -j20^. (8) 

In order to measure gauge- variant observables, each field configuration in the ^/-ensemble 
is placed in a particular gauge, i.e., a gauge tranformation is applied to each configuration in 
the tZ-ensemble which moves the configuration along the gauge orbit to a gauge-equivalent 
configuration satisfying a lattice gauge fixing condition. The simplest gauge choice is defined 
by maximizing 

Y^Tt [f/^ (x) + f/; (x)] (9) 

for each configuration over the class of all gauge transformations. The sum is over all the 
links of the lattice. Any local extremum of this functional satisfies a lattice form of the 
Landau gauge condition: 

5^ \A^ (x + /i) - A^ {x)\ = (10) 

where A^ (x) is a lattice approximation to the continuum gauge field, given by 



U,ix)-UUx) 1 



4 U) = "^ '. ^ - -Tr 



U, (x) - U+ (x) 
2i 



(11) 



A lattice form of Coulomb gauge can be obtained by restricting the sum over fi in Eq. ()1U|) 
above to the spatial dimensions. Other gauge-fixing conditions may also be used |l6|, and 
improvement can be applied to the definition of A^ as well. The global maximization needed 
is often implemented as a local iterative maximization. The issue of Gribov copies arises 
in lattice gauge fixing because such a local algorithm tends to find local maxima of the 
gauge-fixing functional. There are variations on the basic algorithm that ensure a unique 
choice from among local maxima [16[ . 



For analytical purposes, it is necessary to generalize this procedure jll|, so that a given 
single configuration of gauge fields will be associated with an ensemble of configurations of 
(7-fields. We will generate this ensemble using 

S,f [U. g] = Y. ^^^ b"" (^ + ^') U^ (^) 9 (^) + 9^ (^) K {x) g {x + /x)] (12) 

as a weight function to select an ensemble of g fields. The normal gauge-fixing procedure 
is formally regained in the limit a ^ 00. Computationally, this generalized gauge- fixing 
procedure can be implemented as a Monte Carlo simulation inside a Monte Carlo simulation. 
Note that the g fields should be thought of as quenched variables, since they do not effect 
the [/-ensemble. 

The expectation value of an observable O, gauge- invariant or not, is now given by 

(O) = ^ I [dU] e^^t^l^^ / [dg] e^^^'^'^1 O , (13) 

where 

^./[f^] = /Ne"^^'^'< (14) 

Formally, g is just a quenched scalar field. It has two independent symmetry groups, 
Giocai ® G global, so that it appears to be in the adjoint representation of the gauge group, 
but the left and right symmetries are distinct. The generating functional Zgf[U] is a lattice 
analog of the inverse of the Fadeev-Popov determinant [l2| , but there are some important 
differences. Note immediately that Zgf[U] depends on the gauge-fixing parameter a. More 
fundamentally, with the continuum Fadeev-Popov determinant there is the vexing question 
of Gribov copies: what should be done about field configurations on the same gauge or- 
bit satisfying the same gauge condition? The lattice formalism avoids this question. By 
construction, gauge-invariant observables are evaluated by integrating over all configura- 
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tions. Gauge-variant quantities receive weighted contributions from Gribov copies. Thus 
the connection between lattice gauge fixing and gauge fixing in the continuum is not simple. 
There is an apparent confiict between lattice gauge fixing and Elitzur's theorem J18l |. 
which tells us that a lattice gauge symmetry cannot be broken, neither spontaneously nor 
by an explicit symmetry-breaking term in the action. However, gauge invariance implies 
that the lattice gluon propagator is identically zero. Lattice gauge fixing avoids Elitzur's 
theorem by creating a new global symmetry on top of the underlying gauge symmetry. The 
new global symmetry introduced by the lattice gauge-fixing procedure is used to construct a 
proxy for the gauge field with many of the same properties, but transforming as the adjoint 
of the global symmetry rather than the local one. 

III. FORMALISM FOR QUENCHING 

In this section we develop a general formalism for a set of quenched fields, collectively 
designated 02 and a set of unquenched fields, collectively designated 0i. These fields can be 
scalars, spin-1/2 fermions, et cetera. An observable O is coupled to an external source K. 
Using a compact functional notation, the generating functional for O is 

-S2[<t>2;<t>i]+! KO 

We define the generator of 02-connected subgraphs W [K; 0i] by 

so that 

Z[K]= I [rf0i] e-'^['t>.]e'^[K-M-w[^-M ^ (17) 

where the term W [0; 0i] comes from the integration over 02. The difference between the 
quenched and unquenched models is that the quenched model has an additional factor of 
exp (— ly [0; 0i]). This factor cancels all the internal 02 loops, which are absent in the 
quenched model. The two-point function {00) is given by 



Z[K]= I [rf0i] e-^^I'^^l ^t^'^^^L"""^'"^; \ . (15) 



g ■52[</'2;0iJ 



{00) = / [rf0i]e-^^[^i] 



(5^W\ f5W^ ^ 



6K 



(18) 



K=0 

Note that there are two contributions to {00). The first, from 6'^W/SK'^, represents all 
the graphs which are 02 -connected. A graph is 02-connected if it remains connected after 
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(b) 



FIG. 1: Diagrams contributing to {OO). (a) represents the sum of all (/)2-connected graphs con- 
tributing to the two-point function, (b) is an example of a i?i>2-disconnected graph, where the dashed 
line is a (pi propagator. 

cutting any number of 0i lines. The second contribution, from (SW/SK) , represents 02- 
disconnected graphs. See Fig. [T] for examples. 

We can apply the general prescription for quenching to gauge theories. It is important to 
realize that this analysis is aimed solely at understanding the gauge-fixing procedure, and 
does not address the underlying mechanisms for mass generation. For simplicity, we use the 
language of continuum field theories, but the general form of the generating functional also 
holds on the lattice. Our auxiliary fields will be g{x) and its partner g{x) which take on 
values in the gauge group G. The Euclidean action for g is simply 

S, = jF'Tr{D,gyiD^,g), (19) 

where the covariant derivative is Dfj_g = d^g + iA^g. The constant F"^ can be identified 
with a/2N in the lattice theory. The symmetry group is Giocai ® Ggiobau where the gauge 
symmetry acts on the left, and the global symmetry acts on the right. The proxies for the 



gauge fields are the conserved currents associated with the global symmetry: 



j; = -Tr[gT^D,g)^-iDg)T^g+] 



Tr 



T^g+A^g - -T^g^ d ,g 



(20) 
(21) 



If g is expanded around the identity, we see the natural identification of J^ as a gauge-fixed 
form of A"; . 

If we couple a source K'j^ to the currents, then the generating functional is 



Z[K]= / [dA] e 



,,J[ci^]e-W^^ 



n<m 



e ^9 



(22) 



and as before, we define the generator of (y^-connected graphs to be 



^W[K;A] 



The action Sg can be written in the form 



[dg] e 



^Sg+JKJ 



(23) 



Sg = I F'Tr [{d,gy {d,g) + j,A^ + Aj] 



(24) 



where we have defined the gauge-variant currents 



Jm = i {d,^g^) g - ig^d^g . 



(25) 



It is important to note the distinction between J^, a set of currents transforming non-trivially 
under the global symmetry, and the gauge- variant currents j^, which transform under the 
local symmetry, but are invariant under the global symmetry. The generator of (^-connected 
graphs is 



^W[K;A] 



j [dg] exp - j F'Tr [{d,g)^ {d,g) + j.A^ + Aj - K,J^] 



(26) 



Now we have the expression for the proxy of the gauge field propagator 



{j;jl) 



1 



[dA]e 



-Sa 






6W\ f 5W 






where 



5W 



J [dg] J^ exp [- / F^Tt [{d^g^ {d,g) + j,A,] ] 
J [dg] exp [- / F^Tr [id,g)+ (d.g) + j,A^] ] 



(27) 



(28) 
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gives rise to the graphically ^f-disconnected graphs. Similarly, the (^-connected graphs are 
obtained from 



5^W \ _ f [dg] r^Jl exp [- / F^Tr [{d,gY {d,g) + j,A,] ] 
6K-^6Kt) J [dg]exp [- J F^Tr [{d,g)^ {d,g) +j,A,]] 



(29) 



so the total propagator consists of two terms. Note that the apparent mass term for A^ has 
been cancelled out. 

In the case of QED, we can exactly solve for (J^Ju)- Writing g as exp(z^), we have 
Sg = J F"^ {dfj,9 + Afj_) and the current J^ is given by J^ = A^ + d^O. The generating 



functional is 



^W[K;A\ 



j[de]exp - j ^{d,e + A^f + j K,{d^e + A, 



(30) 



where v"^ = 2F^. Then 



^W[K;A] 



exp 



X exp 



—K 



K^A^ 



det {-v^d'') 



-1/2 



[K, 



''A.) i^ (^' 



y2 Q2 



V 



'a:) 



The momentum space form of the propagator is 

fc,,fc 



(J,(A:)J.(-A:)) 



5/.p- 



^ii<^p 



(A,(A;)A,(-A;)) U, 



fcfjfcjy 



i fV/,fC 



H'^v 



(31) 



(32) 



k"^ I \ k"^ I v"^ k"^ 

The first term is the photon propagator, projected onto the transverse subspace, and the 
second term represents a direct contribution from the 9 field. Note that as v goes to infinity, 
the propagator becomes purely transverse. 

It is amusing to note that the above formula is also valid in the unquenched case, where 
no 6 field is introduced. The essential change is that the gauge field acquires a mass v via 
the simplest example of the Higgs mechanism. The propagator for 8^9 + A^ has exactly the 
same form as in the quenched case, but the mass term now causes the A^ two-point function 
to have the explicit form 

{A,{k)K{-k)) = ^-f^^ (33) 

and one easily checks that ( J^ (A;) J^, {—k)) = {A,, {k) A^, {—k)) in the unquenched case. 

Lattice simulations of the ^(1) propagator [8], which correspond to the limit v — > cxd, 
give a non-zero asymptotic mass in the strong-coupling, confining region, and a zero mass 
in the weak-coupling, free field region. In neither region are violations of spectral positivity 
observed, consistent with our results here. 
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IV. SOLUBLE EXAMPLES: MIXING MODELS 

In our previous work [lj|, we analyzed the simplest model displaying spectral positivity 
violation due to quenching: a model of two free, real scalar fields with a non-diagonal mass 
matrix, with one of the fields quenched. Here we generalize this to include the effects of 
mixing in the kinetic terms as well. The Lagrangrian for the scalar mixing model is 

L = \ [{dct^if + m?0?] + i [(902)' + m^0^] - ^^''M2 - e (900 (902) (34) 

where we take 02 to be quenched. If 02 were not quenched, this Lagrangian could be 
rewritten in terms of two free massive fields after a field redefinition. The quenched model 
can be solved by the functional method described above. We have done this in the case 
of . ^ ,n Re. Q. It suffice. Ke.e to note tK. tKe D..o„ .ene. to. t.e ,. p.pa.ato. 
truncated by quenching, is given by 

-^ 9 + -^ 9 {^^^ + ^P^) —r, 9 [l^ + ^P') —0 9 • (35) 

The 02 propagator can also be written in the form 

A ^ B ^ C 

p^ + ml {p^ + mlf p^ + ml 

where A, B, and C are functions of the parameters of the Lagrangian. 

Spectral positivity is violated, because a double pole is a limiting case of a negative-metric 

contribution: 

1 1 r 1 11 

9 I 9 9 I 2 ■ V '/ 



lim 



(p2 + mlf m^m2 ml - m? 
The form in coordinate space is very interesting. In any number of dimensions, we can 
consider propagators using wall sources, i.e., sources of co-dimension 1. This sets the mo- 
mentum equal to zero in all the directions except the direction perpendicular to the wall. 
The coordinate space form of the two-point function is 

2m2 Ami 2mi ^ ^ 

It is the second term, arising from the double pole, which violates spectral positivity. 

Propagators are more complicated when the fields couple to both quenched and un- 
quenched states. Consider a model of the gauge-fixed vector field as a mixture of two 
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FIG. 2: Diagrams that contribute to the propagator in the vector mass- mixing model. Wavy Unas 
represent the Af^ propagator, and straight Unes the Bf^ propagator. 

massive vector fields: i?" which will be quenched, and A" which is unquenched. The 
Lagrangian density can be taken to be 



L = - {d,Al - d^Al) ' + - {d,B: - 9.5^ ' + -m\ (Ay ' + -m\ (5;) ' + /xM^5^ (39) 



with the sum over indices implicit. If both fields were unquenched, then the mass eigenstates 
would be 

— ^^-^ — . (40) 

Consider a field which is a linear combination of A and B, given by kA + XB. Only the 
relative sign of k, and A is important. The use of wall sources greatly simplifies what would be 
a complicated tensor structure, and the relevant Dyson series for the transverse components 
of the propagator is 

K^ A^ ^ kAu^ A^u*^ ,,,. 

I 2 - I - (41) 

p^ + mj p^ + ml {p^ + ml){p^ + ml) {p^ + ml) {p^ + mlf ^ ' 

corresponding to the diagrams shown in Fig. |21 

The propagator can be rearranged in several forms. The factored form 

^,(,_ VA)^(l- VA) (42, 

p^ + rrig \ (p^ + mj) ) p^ ^ m\ \ [p^ + rrig) / 

is helpful because it shows the beginning of the infinite Dyson series in the case where 
both fields are unquenched. The partial fraction form explicitly presents all the information 
available: there is a single pole in p"^ at — mf , a single and a double pole at p^ = —m,^, and 
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three coefficients giving the residues at the poles. 



+ 



mf — 777-2 / P + ^^1 

A/i 



,2 N 2 



mj — 77I2 



+ k] +X^ + k 



2 I ,,2 



p^ + m2 



^^"^ ' (43) 



(mf - m?) (p2 + m2)2 ' 

Note that the residue of the 7771 pole must be positive, but this need not be true for the 
other residues, which come from the quenched fields. 

V. PHASE DIAGRAM OF THE LATTICE MODEL 

The unquenched form of the lattice model is a Higgs model, of a type first analyzed by 
Fradkin and Shenker [l9|. It is very useful to recall their analysis. For a and P small, 
there is a convergent strong-coupling expansion associated with a phase in which the g fields 
are confined into bound states. For a and j3 large, perturbation theory indicates that the 
Higgs mechanism takes place in its most complete form, with no remaining scalar fields. 
At tree level, the only particles are massive vector particles. Naively, there appears to be 
two distinct phases, a confined phase and a Higgs phase. However, the field g is in the 
fundamental representation of the gauge group, and breaks the Z{N) symmetry associated 
with confinement in the pure gauge theory, in a manner similar to quarks. As Fradkin and 
Shenker showed, there is no absolute distinction between the confined and Higgs phases in 
this case, and the two regions of the lattice phase diagram are connected. 

We have studied the phase diagram of the quenched model using lattice simulations, and 
the results are similar to those for the unquenched model. As shown schematically in Fig. |3l 
there is a critical line coming out oi (3 = 00. The origin of the line is the same in both the 
quenched and unquenched models: at /? = cx), the link fields become gauge transforms of 
the identity, and can be absorbed into the g fields. Thus setting j3 = 00 gives a spin model 
with a phase transition between a disordered phase for small a and an ordered phase for 
large a. As (3 decreases, the disordering effect of the gauge fields increases, requiring larger 
a for the phase transition. In both the quenched and unquenched cases, the critical line 
does not completely separate the two phases, but terminates in a critical end point in the 
a — P plane. This line appears to be first order in the quenched case, as is the case in the 
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a 



FIG. 3: Phase diagram in the a — P plane for the SU{2) gauge theory. 



unquenched model. Figure |3] shows the specific heat for the quenched model as a is varied 
for various values of [3 in the case of SU{2). 

As first discussed by Fradkin and Shenker in the unquenched case, the interpretation 
of the massive vector particles of the model depends on which part of the phase diagram 
is being considered. In the confined phase, the vector masses are associated with vector 
bound states of the confined scalars. This interpretation is the same in both the quenched 
and unquenched models. For the unquenched model, the vector particles are interpreted as 
fundamental particles, made massive via the Higgs mechanism. In the quenched model, this 
interpretation is problematic. As we have seen in our continuum treatment of gauge fixing 
as quenching, the term responsible for giving the vector fields a mass cancels out, and the 
gauge fields do not acquire a mass via the Higgs mechanism. The gauge fields can acquire 
a mass via their own self-interactions, however. 
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FIG. 4: Plot of the specific heat for the quenched model as a function of a. 



A. Interpretation as confined theory 



We will interpret the vector multiplet states in the confined region in a manner similar to 
the work of Bardeen et al. on the a^ in quenched QCD J20|, 



21[. The 1]' propagator violates 



spectral positivity in quenched QCD. Within chiral perturbation theory for full QCD, the t]' 
becomes heavier than the particles in the pseudoscalar multiplet by the summation of a mass 
insertion term arising from the anomaly. In quenched QCD, the summation of the Dyson 
series is truncated, giving rise to a double pole in the rj' propagator. The effects of this show 
up in the Oq propagator via loops containing an t]'. The left-hand side of Fig. El^a) shows the 
contribution of a single oq particle to the Oq propagator, and the right-hand side shows an 
equivalent diagram in terms of quarks. Figures El^b) and (c) show loop contributions from 
states containing an t]' as part of the intermediate state. FigureEfb) differs from Fig.^l^c) by 
a single mass insertion. In full QCD, these are the first two terms in a geometric series which 
can be easily summed to give a heavy rj'. In quenched QCD, Fig. ^c) does not contribute, 
because it has an internal quark loop. As shown by Bardeen et al., the bubble in Fig. Efb) 
leads to spectral positivity violation in the oq propagator. 
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FIG. 5: Diagrams contributing to the oq propagator in quenched QCD or the vector propagator in 
the confined phase of a gauge-fixing model. Diagram (c) only contributes in the unquenched case. 

We can take over this argument to case of gauge fixing, where the gauge-fixed vector op- 
erator must be understood as creating bound states in the confined region. In this extended 
analogy, there is a g^Ag vector bound state, whose propagation is represented by Fig. Efa). 
There is also a contribution from intermediate states containing a isoscalar scalar state, as 
in Fig. Efb). The isoscalar scalar thus plays a role in the confining region of the quenched 
theory similar to that of the r]' in quenched QCD. 

Let B{p) stand for the bubble diagram and let the vector propagator be given by V{p). 
We parametrize the coupling between V and B as — /x^. Suppose we look at the propagator 
for an operator that couples with strength n to the bubble and strength A to the vector. 
Then the propagator can be written in the form 



AV(p) 



1 



A/i^ 



V{p) 



K^B{p) 



1 



A/i^ 



K 



V{p) 



(44) 



K - ■ ) l + iJ,^B{p)V{p) 

We have chosen our parameters, including the sign of /i^, to facilitate comparison with 
Eq. (021) • If we identify the pole in the vector propagator V{p) with —rri^, and the pole in 
the resummed bubble n'^B{p)/ (1 + ^^B{p)V{p)) with — -m^, then on shell we have exactly 
reproduced Eq. (J42|) . 
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B. Interpretation in Higgs region 



In order to discuss the region where a is large, we return to the continuum model of 
Sec. IIIH and write g [x) as exp {iH{x)/ F) where H is Hermitian and traceless. We assume 
that for F sufficiently large, we are in the Higgs phase, where g {x) can be expanded around 
the identity. We are only treating g (x) perturbatively, and non-perturbative phenomena in 
the gauge-field sector, such as dynamical mass generation, are not treated here. 

In the unquenched theory, the Higgs mechanism would occur via the term in Sg 



j F'Tr{A,g)^{A,g) = j F'TvAl, 



(45) 



but this term is explicitly canceled in the quenched theory. Internal loops associated with H 
are also cancelled out in the quenched theory, as discussed in Sec. IIIII The gauge- invariant 
current J^, which is a proxy for the gauge field, is 



ja 



Tr 
Tr 



T^g+A,g - -T^g+ d ,g 



T'^[A,~-[H,A^ 



p ^ 



2F2 



[HAH, A, 



2F2 



[H,d,H]\+0{H') 



(46) 



so the current J!J has a term with A't alone. On the other hand, the current which couples 
perturbatively to A^ is 






(47) 



so that Jn and j^ have some common terms. 
The propagator {J^Jl) is given by Eq. (jTTj) 

5^W 



I ja jb\ 



[dA] e 



^Sa 



-) 



6W 



5K15KIJ \5KlJ \5Kl 



5W 



K=0 



The first term represents (^-connected graphs, as shown in Fig. ^a). The second term 
represents all other contributions, including intermediate states involving one or more A^ 
propagators. If we assume that the intermediate state with one A^ propagator dominates, 
then we need only consider the graphs shown in Fig. intb)-(e). The sum of graphs has the 
same structure seen in the confined phase. 
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(b) 



(c) 





(d) 



(e) 



FIG. 6: Schematic diagrams contributing to the vector propagator in the Higgs phase, (a) consists 
of all g-connected contributions, while (b)-(e) are the ^-disconnected pieces with single gluon 
exchange. 

Comparison of Fig. El with Eq. (jl^ immediately suggests that the double pole associated 
with the mass m2 can be understood as arising from Fig.lH^e). This does not fix the origin of 
the mi pole. One possibility is that mi is associated with Fig. IHl^b), in which case mi would 
be a candidate for what we mean by the term gluon mass. As we will see in the next section, 
lattice simulations indicate that the mass mi is too heavy to represent a phenomenologically 
realistic gluon mass. Another possibility is that both the mi and 1712 poles occur in Fig.lHfa), 
as happens in the confined region. From this point of view, the diagrams in Fig.^l^a) and (b), 
which are both (^-connected, could move smoothly into Fig. IHfe). Thus the interpretation of 
the masses is ambiguous in this phase. 

VI. LATTICE RESULTS FOR PROPAGATORS 



We have performed simulations with the action in Eq. Q for the case of SU{2) 20^ 
lattices at /5 = 2.4 and a = {1, 1.2, 1.4, 1.6, 1.75, 2.0, 2.5, 3.0}. The critical line in the (a, (3) 
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plane lies near a = 0.8, so most of these points are well inside the Higgs region of the 
phase diagram. For each value of a, we have performed four independent simulations. Fits 
were carried out for each dataset, and finally combined to give final results for the fitting 
parameters. We begin the analysis in coordinate space to extract the masses from the form 

G{x) = {a + 6m2x)e-"2^ + ce-"'^^ + {x -^ L - x) . (48) 

In order to extract the masses, we first cut the data at x = xq, and fit the term (a + 
bm2x)e~'"^'^^ + {x -^ L — x) to all points x > xq. The value of xq is chosen to achieve the 
best fit, and depends on the dataset, but in all cases Xq > l/rrii. This then gives us a value 
for the mass m2 (and a and b). Once this has been fit, we subtract this from the original 
data set and then fit the remainder to ce"*"^^ + {x —^ L — x) to extract the mass rrii, as well 
as c. 

The coordinate space form of the propagator is sufficient for the determination of the 
two mass parameters. The coefficients k. A, and /i^ are complicated functions of a, b and c 
in the coordinate space propagator. Thus it is easiest to extract the mass parameters mi 
and 1712 in coordinate space, but determine the other parameters in momentum space. After 
Fourier transforming the propagators, we multiply by the factor (p^ + m\Y{p'^ + in\), where 
p^ should be understood as the one- dimensional lattice momentum squared, sin^(p). We 
then fit {p^ + m\Y{p'^ + m\)G{p) to a quadratic polynomial in p^. The coefficients of the 
polynomical can be extracted cleanly, and from these we can determine the parameters from 
Eq. (1411) . yU^, K and A. Empirically, we find that these parameters are only weakly dependent 
on the mass parameters. Finally we have taken the values found with these fitting routines 
and have averaged over the four sets of data for each value of a, where the error bars are 
determined from calculating the standard error as these data are statistically inde pen dent. 

In Figs, m and |H1 we compare our 20^ propagator results with our earlier work |lj] on a 
12^ X 16 lattice at /? = 2.4 for the values a = 1.2 and 2.0 respectively. Errors are smaller than 
the size of the data points. The consistency of the data and the absence of any finite-size 
effects is clear, as expected from the magnitude of the masses. 

In Figs. l9l and 1101 we plot mi and mL2 as a function of a on both 12'^ x 16 and 20^ lattices. 
We note that m.2 appears to be essentially independent of a over the range studied, and 
the results from the two lattice sizes are quite consistent. The mass parameter mi shows 
definite a dependence for smaller values of a. The disagreement in the results for m-i for 
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the two lattice sizes is not statistically significant, but mi is clearly not as well- determined 
as m2- This is to be expected, since mi is much larger than m,2, and is extracted from the 
same propagator. 

The parameters k. A, and fi^ are shown as functions of a in Figs. llH|T!?l The parameter k 
represents the coupling of the lattice gluon operator to the unquenched particle, which has 
mass m,i. It is monotonically decreasing as a function of a, but it is difficult to determine 
if it asymptotes to a non-zero value for large a. The parameter A is the coupling of the 
gluon operator to the quenched particle, which has mass m,2. This parameter is not as well- 
determined as K, but is much larger for large a. The parameter /i^, which represents the 
mixing of the quenched and unquenched particles, show a downward trend with increasing 
a, but with large statistical errors. 

The behavior of the two masses mi and rn.2 as functions of a are completely different 
from the behavior of the vector particle in the unquenched model. As shown in Fig. ^5 
the vector mass m^y in the unquenched model rises as a^'^ over this range of a, all at 
f3 = 2.6. The unquenched propagator fits very well to a single pole form in this region. 
This is completely consistent with the Higgs interpretation of the unquenched model, where 
my is naively proportional to a in physical gauge. Suppose for the moment that we were 
simulating the mixing model instead of an SU{2) gauge theory. Then the parameters mi, 
m,2, and /i from the quenched simulation could be used to determine the mass eigenstates via 
Eq. fjl0|l . This procedure is dubious for the gauge theory, because only one multiplet of vector 
mass eigenstates occurs in the unquenched model, not two. Furthermore, the parameters of 
the quenched and unquenched models are different due to quenching. Nevertheless, if the 
parameters of the unquenched model are used in Eq. (jIUj) . the larger of the two masses is 
commensurate with the single mass observed in the unquenched model. 

Our analysis of the gluon propagator reveals three parameters with dimensions of mass: 
mi, 7712, and ^. As we have discussed above, there are arguments for both mi and m,2 to 
be interpreted as a gluon mass parameter. The mass mi is the mass of the unquenched 
state, and it is natural to assume that it is the gluon mass. However, the mass of the 
lightest glueball in SU{2) is given by md^fo = 3.74 ±0.12 [23, and ay/a = 0.1326, so 
m^cCL = 0.4959. Thus m,i is about 2mG, and is an unlikely candidate for a gluon mass. 

On the other hand, m,2 is slightly larger than uig/^. If we view the lightest scalar glueball 
as being composed of two consituent gluons, then it is natural to identify 1712 as the gluon 
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FIG. 7: The momentum space gluon propagator on the 12^ x 16 and 20^ lattices at a = 1.2. 

mass. Furthermore, m2 is the hghtest mass in the vector channel, and it is logical to identify 
this with the gluon. In the confined phase interpretation of the unquenched model, the 
lightest vector is a bound state of a gluon, scalar, and anti-scalar. In the Higgs phase, the 
lightest vector is simply the massive vector particle formed by the Higgs mechanism. We 
expect that any meaningful gluon mass would be independent of a, and m2 appears to be 
approximately constant over the range of a studied. 

VII. DISCUSSION 



A large number of possible forms for the gluon propagator, in various gauges, have been 



proposed. A brief review is available in Ref. 



We confine ourselves here to a few remarks 



on some of the main themes that have been considered in the literature. We begin by 
noting that there is no evidence from lattice simulations for massless excitations coupled 
to the gluon propagator. In particular, a l/p"^ propagator, which would lead to a confining 
potential from one-gluon exchange [23|, is firmly ruled out. 

The possibility that the gluon propagator might vanish at zero momentum was first 
discussed by Gribov, based on his work on gauge copies j2J]. He proposed that the gluon 
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FIG. 9: Values of the two mass parameters in as a function of a on the 12^ x 16 lattice. The light 
mass, 777-2 is approximately constant as a function of a, while the heavy mass initially decreases 
with increasing a, reaching a somewhat constant value. 
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FIG. 10: Same as Fig. but for the 20^ lattice, 
propagator has complex poles, leading to a gluon propagator of the form 



p^ + m 



4" 



(49) 



This leads to the vanishing of the gluon propagator at p = 0, and a peak in the gluon 
propagator at a non-zero momentum. Others have also explored this possibility and possible 



generalizations 



m 



Zwanziger has proven [27] that the gluon propagator in lattice 
Landau gauge {a = oo) is zero at p^ = in the infinite volume limit. His proof does not 
require the continuum limit, but does require that only configurations in a restricted region 
of configuration space, such as the fundamental modular region, be integrated over. There 
is evidence from some lattice simulations for a peak in the gluon propagator at very low 
wave number |25l . l2q . However, extremely large lattices are required, and the propagator 
does not appear to be trending to zero at p^ = as the volume grows. 

In our lattice simulations, the gluon propagator does not vanish at p^ = and our best fit 
to the data does not indicate a peak to non-zero momentum. However, within the framework 
of our mixing model, it is possible for a peak to occur, depending on the values of various 
parameters. If the gluon propagator is initially rising with p^, it must eventually fall. Thus 
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FIG. 11: K as a function of a. 
the condition for a peak to occur is 

2m\K\^^m\ — 2m\\^ii^ — mlX^m\ + 2m2/tA/i^ — m\\^ii^ — m^K^ > 0. 



(50) 



Our measured parameters are such that this inequaht y n ever holds. We have also fit our 
data with the form used by Cucchieri et al. in Ref. J2^, a generalization of the Gribov 
form. Although a reasonable fit to the data is obtained, the x^ is 50-100 times larger than 
that obtained with our quenched mixing form, and a peak is not always seen. When a 
peak does appear when fitting to the generalized Gribov form, it occurs at very small non- 
zero momenta, directly accessible only on very large lattices. We conclude that a peak is 
a possible feature of the gluon propagator, but the existence of a peak is sensitive to the 
precise form used in fitting. Such a peak may lack fundamental significance, particularly 
since it may only be directly observable on lattices much larger than the scale on which 
confinement appears. 

Our results for the gluon propagator are similar to those of Leinweber et al. [2^ , who fit 
SU{?>) lattice Landau gauge propagators to a variety of possible forms. They obtain a best 
fit with their model A, for which the propagator has the form 



D{P') 



AM^""^ 
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(51) 




FIG. 12: A as a function of a. 



where L is a logarithmic factor mimicking one-loop corrections. Their preferred value for 
their parameter a/, is 2.17^}g. This form is thus rather similar to our quenched mixing form, 
since for a^ = 2, it has a triple pole. Of course, we have two mass scales, m\ and m2, playing 
roles in the denominators. 

We have not explored the interesting possibility that the gluon propagator has an anoma- 
lous dimension, first suggested by Marenzoni et al. J30|. This behavior could occur in our 
formulation. For example, the bubble diagram discussed in Sec. |V] could give this behavior. 
However, high precision data as well as a detailed functional form for fitting would be needed 
to test for this behavior. 



VIII. CONCLUSIONS 



We have shown in detail how the quenched character of lattice gauge fixing leads to 
violations of spectral positivity in non-Abelian gauge theories. Although we have focused 
on lattice Landau gauge, this problem is likely to occur in any lattice gauge fixing scheme in 
which the gauge fixing variables are determined subsequent to the generation of lattice gauge 
field configurations. A key step in our analysis has been the generalization of conventional 
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lattice gauge fixing, where one configuration is chosen on the gauge orbit, to a forniahsm 
including a gauge-fixing parameter, similar to continuum gauge fixing. This parameter 
controls the weighting of configurations along the gauge orbit, and conventional lattice 
gauge fixing is formally recovered as a limiting case. We have not examined lattice gauge 
fixing schemes in which the gauge field integration is restricted to a subset of the entire 
configuration space, e.g., the fundamental modular domain, but the relation of these to 
standard computational practice is unclear. There are some lattice gauge fixing schemes, 
such as lattice Laplacian gauge J3l| which are defined by a complicated algorithmic process, 
and may not generalize in a manner similar to lattice Landau gauge. Nevertheless, we believe 
that the quenched character of these algorithms is the root of spectral positivity violation 
there as well. 

We have focused here on the form of the gluon propagator, the extraction of masses, their 
interpretation and dependence on the gauge- fixing parameter a. The issue of the continuum 
limit for gluon propagator masses has not yet been resolved. In the gauge-invariant sector, 
physical masses scale as the gauge coupling (3 is taken to infinity in the manner prescribed 
by the renormalization group. In principle, it is also necessary to adjust a to obtain a 
continuum limit of gauge- variant quantities. In other words, the running of the gauge-fixing 
parameter must be taken into account. This is similar to the running of the gauge- fixing 
parameter in covariant gauges in the continuum. If a meaningful gluon mass, independent 
of a for large a, can be extracted, it might scale correctly in the continuum limit without 
tuning a. As we have observed, the mass 1112 is approximately independent of a for large a 
aX (3 = 2.6, and has a value phenomenologically consistent with a constituent gluon mass. 
Further study is required to see if m2 and mi survive in the continuum limit. Because the 
standard gluon operator creates a mixture of states, it is likely that a variational analysis 
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mi. 

We have used the group SU{2) to study gluon properties at zero temperature. Perhaps 
the most important application of lattice gauge fixing is the measurement of electric and 
magnetic masses in the deconfined phase. We are extending our work to finite temperature 
SU{2), generalizing the work of Refs. [3|, y, |5| to the case of finite a. However, it is SU{^) 
which is of principal interest |6|. Because gauge fixing can be carried out any time after lattice 
field configurations are generated, it is relatively easy to study gauge fixing on available 
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unquenched SU{3) configurations. In that case, it would be possible to also examine the 
dependence of the quark mass on the gauge parameter a. 
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